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This paper gives an overview of recently developed model for the QCD analytic invariant 
charge. Its underlying idea is to bring the analyticity condition, which follows from the 
general principles of local Quantum Field Theory, in perturbative approach to renor- 
malization group (RG) method. The concrete realization of the latter consists in explicit 
imposition of analyticity requirement on the perturbative expansion of function for 
the strong running coupling, with subsequent solution of the corresponding RG equa- 
tion. In turn, this allows one to avoid the known difficulties originated in perturbative 
approximation of the RG functions. Ultimately, the proposed approach results in quali- 
tatively new properties of the QCD invariant charge. The latter enables one to describe 
a wide range of the strong interaction processes both of perturbative and intrinsically 
nonperturbative nature. 
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ning coupling. 
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1. Introduction 

Theoretical analysis of strong interaction processes rests on the quantum non- 
Abelian gauge field theory, namely, Quantum Chromodynamics (QCD). This theory 
is dealing with quarks and gluons. The former are the constituents of hadrons, while 
the latter are quanta of a massless gauge vector field, which provides interaction 
between quarks. Gluons and quarks carry a specific quantum number, which plays 
a key role in hadron physics. It was first introduced in the mid 1960s in Ref. T, and 
later it was name "color" . The comprehensive theoretical and experimental in- 
vestigations revealed that the strong interaction possesses two distinctive features. 
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First, the interaction between quarks inside the hadrons weakens with increasing 
the characteristic energy of a process. In other words, the invariant charge*^ vanishes 
at large momenta transferred, representing the so-caUed asymptotic freedom of the 
theory. Second, free quarks and gluons have not been observed, at least until the 
present time. Here we are dealing with the color confinement. 

These two phenomena are related to hadron dynamics at different energy scales. 
The asymptotic freedom takes place at large values of energies, or in the so-called 
ultraviolet (UV) region. It corresponds to small (r < 0.1 fm) quark separations. On 
the contrary, the color confinement is related to the low-energy, or infrared (IR) 
domain, that corresponds to large quark separations, namely, r > 1.0 fm. 

Theoretical description of the asymptotic freedom was first performed in the 
mid 1970s in Refs. O These papers gave rise to wide employment of perturbative 
calculations in Quantum Chromodynamics. Indeed, if the value of invariant charge 
is small, then one can parameterize unknown quantities by perturbative power 
series in the running coupling. On the one hand, this drastically simplifies analysis 
of hadron interactions in the UV region. On the other hand, the hadron dynamics 
at low energies, in particular, the quark confinement, entirely remains beyond the 
scope of perturbation theory. As a rule, description of the QCD experimental data 
in the IR region requires invoking other nonperturbative approaches, for instance, 
phcnomcnolo gical p otential models fRefs. I^12|l. sum rules method (Refs. [T5hlt)|l . 
string models p^ l -*^^ ! bag models,^! variational perturbation theory (Refs. EDHSl) 
and lattice simulations fRefs. Elt429(l . 

The renormalization group (RG) method plays a fundamental role in the frame- 
work of Quantum Field Theory (QFT) and its applications. It was developed almost 
half a century ago in Refs. '301 and 13 II (see also Rcf. 32). Nowadays, theoretical anal- 
ysis of hadron dynamics relies, in main part, upon the RG method. Meanwhile, 
the construction of exact solutions to the renormalization group equation, which 
is a rather nontrivial mathematical problem itself, is still far from being feasible. 
Usually, in order to describe the strong interaction processes in asymptotical UV 
region, one applies the RG method together with perturbative calculations. In this 
case, the a priori unknown renormalization group functions are parameterized by 
the power series in the QCD running coupling. Ultimately, this leads to approxi- 
mate solutions of RG equation, which are commonly used for quantitative analysis 
of experimental data. However, employment of perturbation theory in the frame- 
work of renormalization group method generates unphysical singularities of such 
solutions. The latter are in contradiction with the basic principles of local Quan- 
tum Field Theory. For example, solution to RG equation for invariant charge gains 
unphysical singularities due to perturbative approximation of the /? function. At 
the one-loop level there is the so-called Landau pole, and account of the higher-loop 
corrections just introduces additional singularities of the cut type into expression 

"Sometimes invariant charge a(g^) = 9^(9^)/(47r) is also called the (strong) running coupling 
of QCD. 
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for the running coupling (see [Appendix Aj l. 

An effective way to overcome such difhcuhies was proposed"^-' in the late 1950s 
by P.J. Redmond. It consists in invoking into consideration the analyticity require- 
ment, which follows from the general principles of local QFT. This prescription 
became the underlying idea of the so-called analytic approach to QFT, which was 
formulated in the framework of Quantum Electrodynamics (QED) by N.N. Bogoli- 
ubov, A. A. Logunov, and D.V. Shirkov.34s Really, the QED running coupling is 
proportional to the transverse part of the dressed pho ton propagator. Therefore, 
in accordance with the basic principles of local QFT the spectral representation 
of the Kallen-Lehmann type holds for it. The latter implies the definite analytic 
properties in variable for the invariant charge. Namely, it should be the analytic 
function in the complex g^-plane with the only cut along the negative*^ semiaxis of 
real . 

In the general case the QCD invariant charge is defined as a product of propaga- 
tors and a vertex function (see, e.g., Ref.|SS|). Therefore, one might pose a question 
concerning the analytic properties of this quantity. This matter has been examined 
in Ref. I36L and it was shown therein that in this case the integral representation 
of the Kallen-Lehmann type holds for the running coupling , to o. Proceeding from 
these motivations, the analytic approach was lately extended^ to Quantum Chro- 
modynamics by D.V. Shirkov and I.L. Solovtsov, and applied to the "analytization" 
of perturbative series for the QCD observables'^ (see reviews EH and Refs. l39H42l for 
details). The term "analytization" means here the restoring of the correct analytic 
properties in the q^ variable of a quantity under consideration by making use of the 
Kallen-Lehmann integral representation 



{a(<z^)} 







^rf- (1) 

a + q-^ 



In this equation the spectral function Q[a) is determined by the initial (perturba- 
tive) expression for a quantity in hand 

e(cr) = lim [a(-ct - ie) - A(-(T + ze)j , cr > 0. (2) 

ZTTl £— ►0+ L J 

It is worth mentioning also the appealing features of the analytic approach to 
Quantum Field Theory, namely, the absence of unphysical singularities and fairly 
well higher loop and scheme stability of outcoming results. 

The description of hadron interaction at low energies remains a long-standing 
challenge of Quantum Chromodynamics. Some insight into many basic problems 
of elementary particle physics can be gained form constructing the models, which 
adequately account both the perturbative and intrinsically nonperturbative traits 



'^A metric with signature (—1, 1, 1, 1) is used, so that positive corresponds to a spacehke mo- 
mentum transfer. 

•^This method has been named the analytic perturbation theory (APT) thereafter. 
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of hadron dynamics. Certainly, such studies are useful for further development of 
the strong interaction theory. Besides, the description and interpretation of many 
experimental data can be performed only by invoking such models. Thus, the results 
obtained in this field could expand our scope of knowledge concerning the nature 
of fundamental interactions. 

This paper overviews a new model for the QCD analytic invariant charge. The 
model is formulated in detail with a specific attention to the properties of the 
analytic running coupling and to verification of the model's self-consistency. The 
principally new way of involving the analyticity condition into the RG formalism is a 
distinguishing feature of this model. Namely, the analyticity requirement is imposed 
on the perturbative expansion of the renormalization group (3 function. Eventually 
it leads to a number of profitable traits of this running coupling, that enables one 
to describe a wide range of the strong interaction processes both of perturbative 
and intrinsically nonperturbative nature. It is worth noting also that this model 
for the strong running coupling has no adjustable parameters. Therefore, similarly 
to perturbative approach, Aqcd remains the basic characterizing parameter of the 
theory. 

The layout of the review is as follows. 

In Section[21the model for QCD analytic invariant charge is discussed. The key 
point here is the restoring of the correct analytic properties of perturbative ex- 
pansion for the f3 function, and the subsequent solution of the corresponding RG 
equation. The explicit expression for the one-loop analytic invariant charge (AIC) 
is obtained. It is proved to have qualitatively new properties. In particular, the 
analytic invariant charge has no unphysical singularities, and it incorporates UV 
asymptotic freedom with IR enhancement in a single expression. The AIC is found 
to be in agreement with the results of various nonperturbative studies of hadron dy- 
namics in the infrared domain. For example, this strong running coupling explicitly 
reproduces the recently discovered symmetrjl^ related to the size distribution of 
instantons. The consistency of the considered model with general definition of the 
QCD invariant charge is proved. The analytic running coupling is also represented 
in explicitly renorminvariant form, and the properties of the relevant (3 function 
are studied. A number of symmetry relations for the analytic invariant charge and 
the corresponding /? function, which establish a link between the ultraviolet and 
infrared domains, are derived. The different ways of incorporating the analyticity 
requirement into the RG formalism are discussed. 

Section 13 is devoted to investigation of the analytic invariant charge at the 
higher-loop levels. The integral representation for this running coupling is derived 
herein. Similarly to the one-loop level, analytic invariant charge is shown to have 
no unphysical singularities. The developed model possesses a good higher loop and 
scheme stability. The investigation of the corresponding P function revealed that 
AIC has the universal asymptotics in both the ultraviolet and infrared regions, ir- 
respective of the loop level. Further, the model in hand is extended to the timelike 
domain. The difference between values of the analytic running coupling in respec- 
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tive spacelike and timelike regions is found to be considerable for intermediate and 
low energies. Apparently, this circumstance must be taken into account when one 
handles the experimental data. The result obtained in this section confirms the hy- 
pothesis due to Schwinger concerning the (3 function proportionality to the relevant 
spectral density. 

The applications of the analytic invariant charge to study of a number of strong 
interaction processes are gathered in Section 01 In particular, the static quark- 
antiquark potential constructed by making use of the AIC is proved to be confining 
at large distances. At the same time, at small distances it has the standard behav- 
ior originated in asymptotic freedom. Further, the developed approach is applied 
to description of gluon condensate, inclusive r lepton decay, and electron-positron 
annihilation into hadrons. The congruity of estimated values of the scale parame- 
ter Aqcd testifies that the analytic invariant charge substantially incorporates, in 
a consistent way, both perturbative and intrinsically nonperturbative aspects of 
Quantum Chromodynamics. 

In the Conclusions (Section O the properties of the analytic invariant charge 
are summarized and the basic results are formulated in a compact way. 

The auxiliary materials are collected in the appendixes. In |Appendix A| the per- 
turbative solutions to the RG equation for QCD running coupling are derived at 
different loop levels. In [Appendix B| thc Lambert W function is briefly described. 
The function N(a), which essentially simplifies investigation of the /3 function cor- 
responding to the one-loop AIC, is defined and studied here. [Appendix C] contains 
explicit expressions for the spectral functions used in construction of the analytic 
invariant charge. In [Appendix D[ the coefficients specifying the expansion of the 
quark-antiquark potential are presented. 



2. The QCD Analytic Invariant Charge 
2.1. Formulation of the model 

First of all, let us consider the renormalization group equation for the strong running 
coupfing (see, e.g., R.efs. r^and .441: 



^ln[.9^(M)] 
d In fi^ 



/3(5(m)). 



(3) 



In the framework of perturbative approach, the (3 function on the right-hand side 
of this equation can be represented as a power series 



/3o 



9M. 
167r2 



167r2 



(4) 



where /3o — 11 — 2nf/3, (3i — 102 — 38nf/3, and rif is the number of active quarks. In- 
troducing the standard notations ct^ifJ.'^) = .g^(/x)/(47r) and a(/i^) — a(/i^) /3o/(47r), 
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one can reduce the RG equation Q at the ^-loop level to the form 



d In /i2 '-'3 



3=0 



/3o 



(5) 



It is well-known that the solution to this equation has unphysical singularities at 
any loop level. So, there is the Landau pole at the one-loop level, and account of 
the higher loop corrections just introduces the additional singularities of the cut 
type into expression for the invariant charge.'^ However, the fundamental principles 
of local Quantum Field Theory require the invariant charge a{q^) to have a definite 
analytic properties in the variable. Namely, there must be the only cut along the 
negative semiaxis of real (see Ref. I57|l. 

The solution to the renormalization group equation Q gains unphysical singu- 
larities due to perturbative approximation of the P function Q . Such a representa- 
tion of the right-hand side of the RG equation proves to be justified in the ultraviolet 
region, where the asymptotic freedom takes place. However, the approximation Q 
does not ultimately lead to the correct analytic properties of the QCD invariant 
charge. In fact, since the strong running coupling a{q^) is of a fixed sign, the con- 
dition that invariant charge has the only left cut in q^ means that the right-hand 
side of the RG equation 0, as a function of q^, has no singularities in the region 
of positive q^. Apparently, the perturbative expansion violates this condition. 

In the framework of developed mode^ the analyticity requirement is im- 
posed on the P function perturbative expansion for restoring its correct analytic 
properties. This leads to the following RG equation for the QCD analytic invariant 
charge: 



d In /i^ 




In this equation denotes the ^-loop analytic invariant charge, ai^\^^) is the 

perturbative running coupling at the £-\oop level, and the braces { ■ ■ • mean the 
analytization of the expression within them by making use of the Kallen-Lehmann 
spectral integral (Q) (see also R.efs. .i47;,andi48|l. 



2.2. The one-loop analytic invariant charge 

Since the proposed model rests on invoking the analyticity requirement into the 
renormalization group formalism, one may anticipate that the strong running cou- 
pling acquires qualitatively new properties. Let us proceed to construction of solu- 
tion to the RG equation ((HJ, restricting ourselves to the leading loop approximation 
at this stage. 



Solutions to equation JSJ at different loop levels are considered in detail in [Appendix A| 
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At the one-loop level the renormalization group equation ® takes the fornJSEi 



. (7) 

^ J an 



d In/i^ 

Here ain (m^) is the one-loop analytic invariant charge and ai^\^'^) denotes the 
perturbative running coupling at the one-loop level. Taking into account Eq. 
one can represent the right-hand side of the RG equation lO as follows 

a In /i^ J a + fj,'^ 



where Tlf-^^a) = [ln^((T/A^) -I- 7r^]~^. Upon the integration with respect to a and 
the introduction of the dimensionless variable = p? j b^^ ^ Eq. © acquires the 
form 

dln[aW(^^A^)] 



1 1 

d In In 1 — 1^2 

Further, integrating this equation with respect to in finite terms, one gets 

3in^(g^) _ z - 1 zq Inzo 



iin^go") ^^^^ 



1 ' 



(9) 



(10) 



where z = j ts?' and zq = q^/h?. Thus, the QCD analytic invariant charge at the 
one-loop level has the form" 

Po z mz A^ 

It is worth mentioning also that for this strong running coupling the Kallen- 
Lehmann representation 

oo 

a^Kq^) ^ -T da 12 



holds, where p'"'^^(cr) is the one-loop spectral density 

^(i)(^)=fl+n^^. (13) 



y In^ (7 + TT 

Figure n depicts the analytic invariant charge a'-tniq^) — ain (g^) /3o/(47r) to- 
gether with the one-loop perturbative running coupling ai^\q'^) = 1/ ln((7^/A^). 
The obtained solution to the RG equation pil) possesses a number of profitable 
features. First of all, the integral representation of the Kallen-Lehmann type p2|l 

'^Equation llli has been derived by employing the condition ain{9^) ~* c*s^'{9^)i when — > oo 
(see also discussion in Subsection lij.li . 
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Fig. 1. The one-loop analytic invariant charge ain'(g'^) defined by Eq. 111! (solid curve) and its 
perturbative analog as^\q^) (dashed curve), z = q^/A^. 



implies that the analytic invariant charge ainl^^) has the only cut^ along the 
negative semiaxis of real q"^. Thus, the prescription proposed ultimately leads to 
elimination of the Landau pole in a multiplicative way. One has to note that the 
analytic invariant charge Hll() incorporates the ultraviolet asymptotic freedom with 
the infrared enhancement in a single expression. In particular, this trait will play 
an essential role in applications of the elaborated model, see Section0]further. It is 
worth emphasizing also that the developed model contains no adjustable parame- 
ters, i.e., similarly to perturbative approach, Aqcd remains the basic characterizing 
parameter of the theory. 

A crucial insight into the nonperturbative aspects of the low-energy hadron dy- 
namics can be provided by the relevant lattice simulations. One should mention 
here a recent lattice investigation of the topological structure of the SU(3) vacuum 
by the UKQCD CoUaboration.EHI It was revealed therein that the size distribution 
of instantons D{p) has a conspicuous peak shape, and it is severely suppressed at 
large scales p. In turn, this is compatible neither with perturbative results, nor 
with the freezin g of the strong running coupling to some constant value at large 
distances.ESEOI Remarkably, the analytic invariant charge piB is proved to repro- 
duce explicitly a special kind of symmetry related to this issue^*^ (see Subsection l2.[{l 
for details). Meantime, the strong running coupling itself can also be studied on the 
lattice. In pa rtic ular, a recent investigation of its low-energy behavior by ALPHA 
CollaboratioJSl testify to the infrared enhancement of the QCD invariant charge. 
The study of the Schwinger-Dyson equations is another source of the nonpertur- 
bative information concerning the strong interaction issues. It is worth noting here 

'The function a{z) = (z — l)/(z Inz) has a smooth behavior in the vicinity of the point z = 1, 
namely a(l + e) ~ 1 - e/2 + 0{e^), £ -» 0. 
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that the behavior of the QCD invariant charge, similar to Eq. is found to be 
in agreement with the solution of the Schwinger-Dyson equations (see papers 1521 
and references therein for details). 

In general, the QCD invariant charge is defined as the product of the corre- 
sponding Green functions and vertexes. For example, in the transverse gauge the 
following definition 

a{q^) = a{^?)G{q^)g\q^) (14) 

takes place (see, e.g., Ref. OSJ). In this equation a{q^) is the QCD running coupling, 
Q!(/i^) denotes its value at a normalization point /i^, G{q^) and g{q^) are dimension- 
less gluon and ghost propagators, respectively. It is essential to note here that the 
Green functions possess the correct analytic properties in the q^ variable^ before 
the RG summation. Indeed, in massless theory the Green function can be written 
as the perturbative power series in a(/z^) with the coefficients being a polynomial in 
ln((7^/^^). Obviously, in this case the integral representation of the Kallen-Lehmann 
type must hold for any finite order of perturbative expansion for the Green func- 
tion.'^ However, the RG summation leads to violation of such a property of the 
Green functions (e.g., in the simplest one-loop case the Landau pole appears). As 
it has already been mentioned, the analytic approach to QCD seeks to restore the 
correct analytic properties of the relevant quantities. It turns out that for the con- 
sistency of involving the analyticity requirement with the definition of the invariant 
charge (|14|l . one has to apply the analytization procedure to the logarithmic 
derivative of the Green function d h\G{q^) / d h\q^ . Really, if the spectral function 
of the Kallen-Lehmann representation for the Green function G{q^) is of a fixed 
sign (that is true for the leading orders of perturbation theory), then G{q^) has no 
zeros in the complex g^-plane. Hence, the derivative d lnG{q^)/d \nq^ can also be 
represented as the spectral integral of the Kallen-Lehmann type (cf. equations ((S)) 
and H12(l for the one-loop level, and Eqs. 137|l and 142|l for the higher loop levels). 

In the framework of perturbative approach the Green functions considered above 
have the following form: g{q'^) = (1/ln z)''9 and G{q'^) = (l/ln zY° , where z = 
q^ / , dg and da denote the corresponding anomalous dimensions. For the latter 
the relationship 2dg + da = ^ holds, that plays the key role in definition p4(l . 
Applying the analytization procedure to these functions in the way specified above, 
one arrives at the following result: g^niQ^) = [(^ ^ l)/(-2 Inz)]'*" and G^j^{q'^) = 
[{z — l)/(z Inz)]'''^. Therefore, in this case invoking the analyticity condition does 
not affect' the definition of the invariant charge (|14|l . Thus, we infer that the analytic 
running coupling is consistent with the general definition of the QCD invariant 
charge (see also discussion of this matter in Ref. I54|l . 

^Namely, there is the only cut along the negative semiaxis of real q^. 

''One should note that the relevant spectral function here can not be directly identified with the 
spectral density. 

'It is interes ting to note that the similar situation takes place in the nonperturbative a-expansion 
method alsoESlSl 



10 



2.3. Properties of the one-loop analytic invariant charge 

Let us address now the renorminvariance of the analytic running couphng Hll|l . In 
general, any expression for the invariant charge makes sense only if the relevant 
definition of the parameter Aqcd is provided. Otherwise, the running coupling may 
not be a renorminvariant quantity at all. In fact, the renormalization invariance 
of a solution to the RG equation (0) is ensured by a certain relation between the 
parameter Aqcd and the value of the running coupling a{^'^) at a reference point 
(see Eqs. (jA.6|) and (|A.10|I also). However, the way of introduction of the parameter 
Aqcd, considered in the previous section, blurs such a relation. Nevertheless, in the 
framework of the model in hand, the latter can be recovered by solving the equation-i 



a 



(i; 



(9^ 



(15) 



with respect to parameter A. In [Appendix B| equations of this type are studied in 
detail. Thus, the solution to Eq. H15|l can be explicitly written down in terms of the 
Lambert W function which is defined by the relation'* 



Wkix) exp Wk{x) = X. 
For our purposes it is convenient to introduce the function N{a) 



N{a) 




- exp I 

a \ a 



(16) 



(17) 



In Eqs. Hlt)|) and (|17|l k denotes the branch index of the Lambert W function. It is 
straightforward to verify that the solution to Eq. (|15|) has the form 



47r 



(18) 



Let us note here that the right-hand side of this equation tends to the one-loop 
perturbative form when fx^ ^ 00 (see Eqs. HA.6|I and (|B.16|I also): 

A2=^2g^p . (19) 

[ ^oa(M)J 

Thus, the renorminvariant expression for the analytic running coupling Hll|l is the 
following: 



a 



(1) 



A-K z - \ 



1 

A2 



4^, A^ 



? N 



(3q z \nz' 

where the function N{a) is defined in Eq. H17|) (see also Ref. iSGjfor details). 



(20) 



^Equation I15i is identical to the normalization condition g^{l,g) = for the QCD invariant 
charge g'^ {q^ / fi^ , g) / (Att) = (see, e.g., Ref. l55l . 

''The properties of the Lambert W function 11(51 and the function N{a) 1171 are described in 
detail in [Appendix B| 
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In the framework of the elaborated model the analyticity requirement is im- 
posed on the P function perturbative expansion for restoring its correct analytic 
properties. Apparently, the form of the 13 function in this case is affected by the ana- 
lytization procedure. Thus, it is of a particular interest to study both the /? function 
itself and its properties within the approach developed. 

In general, the (3 function corresponding to the one-loop analytic invariant 
charge is defined by' 

d In a(/i^) 



a 



(1) 



(21) 



d In ' 

The function N{a) (see Eq. (^TJ) enables one to obtain the explicit expression for 
the right-hand side of Eq. (j^Jl: 



1 - N{a)/a 



(22) 



ln[7V(a)] ■ 

In order to examine the asymptotics of the (3 function 122|) it is worth employing the 
Eqs. (|B.16p - ljB.18p . It turns out that for small values of the running coupling a(/i^) 
equation H22I) coincides with the well-known perturbative result 



[1 ( 




-exp 




a V 





0. 



(23) 



The second term in this equation implies the intrinsically nonperturbative nature 
of the (3 function (|22|l . Then, for large values of a(/i^) one has 



-l + O 



ln(ln a) 



(24) 



(lna)2 

Such a behavior of the (3 function leads to the infrared enhancement of the QCD 
invariant charge, namely (up to a logarithmic factor) a(g^) ~ h? / , when — > 0. 
It is worth mentioning also that 



0{e% 



0, 



(25) 



1 e 
2^6 

i.e., the function (|22|l has a smooth behavior at the point of interchange between 



the real branches of the Lambert W function. Moreover, /3in (a) < for any non- 
negative value of aip^), that ensures the asymptotic freedom of the theory. 

Figure El presents the (3 function H22(l together with its perturbative analog 



which corresponds to the one-loop running coupling a^^ {q 



A" {a) 

l/ln(q^/A^). In particular, this plot visualizes the reproduction of perturbative 
limit (|23|) within the model elaborated. 

Thus, the proposed way of involving the analyticity requirement into RG formal- 
ism eventually leads to the qualitatively new features of the (3 function (|22|l . Namely, 
it incorporates the asymptotics related to both the UV asymptotic freedom ^2'A\ 



'In Ref. 1561 the definition /3(a) = da{p?)/d In/x^ was used. 
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Fig. 2. The (3 function I22i (solid curve) corresponding to the one-loop analytic invariant 
charge llli . The relevant perturbativc result /3s^'(a) = —a is shown by the dot-dashed line. 

and the IR enhancement H24() of the QCD invariant charge. In particular, such be- 
havior of the (5 function corresponds to the so-called "V~scheme" (see Refs. 
and l57fl . which is widely used in studies of the quark confinement. 

Let us address now the symmetries of the analytic invariant charge ain (9^) and 
the corresponding (3 function /jin^ (a) , which relate these quantities in the ultraviolet 
and infrared domains. Obviously, such relations are decisive, since they provide us 
with insight into the intrinsically nonperturbative features of the strong interaction. 
The symmetry of this kind for the analytic invariant charge Ulljl was found in 
Refs. island El 

= f (26) 

It is interesting that there are similar symmetry relations for the (3 function (|22(l 
also. In particular,™ 

+/5(5i^nHAV9')) = -1 (27) 

and 

pia^Xq^)) /3(5i'J(AV9')) 

where /3(a(/x^)) = d In a(/i^)/d In //^. In fact, there is a perturbativc analog of 
Eq. (PSl) . namely (3{a^\q^)) /ai'\q^) = — 1, which holds for all non-negative . 

From the general viewpoint, one might anticipate that these symmetries, inter- 
relating the hadron dynamics at short and long distances, could be revealed also 

■"Equation 1271 can be deduced from Eq. 1261 . 
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in some intrinsically nonperturbative strong interaction processes. It is worth men- 
tioning here the recently disc over ed conformal inversion symmetry related to the 
size distribution of instantonsP^ At the classical level, the whole instanton sector 
is proved to be invariant under the four-dimensional conformal transformations!^ 
In particular, the coordinate inversion 

^2 

Xf^ ^ x^ — (^'^) 

transforms the SU(2) Yang-Mills instanton solutiorPSEOIof size p in singular gauge 
to the anti-instanton solution of size Pg/p in regular gauge. However, at the quan- 
tum level it is rather difficult to expect that the instanton sector is still invariant 
under the conformal group transformations since the scale invariance is broken via 
the regularization and renormalization. Nevertheless, there are evidences that the 
invariance under the coordinate inversion H29|l , which has a number of important 
physical implications, exists, in a certain form, at the quantum level, too. Thus, the 
essential result in this field was obtained in Ref . 021 where it was observed that the 
high-quality lattice simulation datc&2^^ for the quantity 



D-\p) = D{p)[^] , D{p)=h 
Po 



2tt 
a{p) 



exp 



27r 
a{p) 



(30) 



appear to satisfy with high precision the conformal inversio n sy mmetry D'~^^{p) — 
D^'{p^/p). Here D{p) is the size distribution of instantons,!^ p is the radius of 
an instanton, a{p) denotes the strong running coupling, po is the peak position 
of D'~'^{p), & is a constant factor depending on the gauge group, and Nc is the number 
of colors. This symmetry of D'^^{p) is evidently originated in the mentioned above 
conformal invariance of the instanton sector at the classical level, and it turns out 
to be at most weakly broken at the quantum level. 

Apparently, the conformal inversion symmetry D'^^(p) — D'^^[p^/ p) imposes a 
certain constraint onto the QCD invariant charge, leading to an important relation 
between the properties of the a{p) at small and large distances. In fact, one can 
derive the expression for the nonperturbative strong running coupling by employing 
the symmetry of this kind. Following thi s w ay, the analytic invariant charge Hll|l 
has recently been rediscovered and provecP^ to reproduce explicitly the conformal 
inversion symmetry of D'^^{p). This remarkable result opens up an alluring op- 
portunity to build up a new general method of consistent description of the both 
perturbative and intrinsically nonperturbative aspects of hadron dynamics towards 
the infrared regime. 



2.4. Other models 

Several decades ago, when the analytic approach to Quantum Electrodynamics 
was developing, it was noteci^ that there is no unique way of incorporating the 
analyticity condition into the RG formalism. There was also an attempll^to resolve 
this ambiguity by involving an additional condition, in particular, by making use 
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of the equations of motion. Following this way the photon propagator has been 
derived^"' the unphysical singularity being removed in a multiplicative way (i.e., 
similarly to Eq. Hll|l l. Furthermore, the presence of nonperturbative terms both 
in the running coupling itself and in the corresponding spectral density 
is another common feature of the model in hand and the approach considered in 
Ref. inSI(see also discussion of this issue in R.efs. 1461 1541 and . 

The identical situation arises in Quantum Chromodynamics, namely, different 
models for the strong running coupling can be proposed in the framew ork of the 
analytic approach to QCD. As it was mentioned above, our modePSESlig based on 
recovering the correct analytic properties of the perturbative expansion of the RG 
/3 function 0. Meanwhile, in the original model due to Shirkov and Solovtsov"''^ the 
analytization procedure 1^ is directly applied to the perturbative invariant charge 
(see also reviews SSj . At the one- loop level it gives 

^2 



An 



1 

In z 



1 



1- z 



Si 
A2- 



(31) 



From the general point of view, one can anticipate the different properties of these 
two running couplings. On the one hand, both these models have no unphysical 
singularities and reproduce perturbative limit in the ultraviolet region. On the 
other hand, there is a difference between them in the infrared domain. Namely, 
the analytic invariant charge possesses the IR enhancement, while the running 
coupling (|31|l freezes to the universal constant value Q!ss(0) = Att/Pq. One has to 
emphasize here that both these models have no adjustable parameters, so they are 
the "minimal" ones in this sense. A brief discussion of this issue can also be found 
in Refs.EHlElandinil 

It is worth mentioning some other models for the strong running coupling within 
the analytic approa ch to QCD. Let us start with the model developed by Alekseev 
and Arbuzov !*^ I'^'^l By making use of a special solution to the Schwinger-Dyson 
equations for the gluon propagator, these authors proposed the following expression 
for the QCD running coupling: 



47r 



1 
In z 



1 



1 -c 



I — z z z + m^/A^ 



A2' 



(32) 



where rUg is the gluon mass and c denotes a dimensionless parameter fixed by the 
phenomenological value of the gluon condensate. Likewise the analytic invariant 
charge (lllll . the running coupling H32|) has enhancement in the infrared domain. 

It is interesting to note here that another similar model for the QCD invariant 
charge 



An 
To 



1 

In z 



1- z 



T_ 

A2 



(33) 



has been put forward in Ref. 65 for the analysis of the lattic e simulation data on 
the low-energy behavior of the QCD running couplingE^Eni The model (|33f) also 
possesses the infrared enhancement. 
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By making use of a certain phenomenological reasoning, Webber suggested the 
strong running coupling of the foUowing forrn"^ 



-,(1) 



(9 ) = 

Po 



1 

In z 



1 z + h fl 



1- z 1 + 5 \z + c 



A2' 



(34) 



with a specific choice of the parameters 6 = 1/4, c = 4, and p = 4. On the contrary, 
this model has infrared finite value ^^•'(O) ~ 2tt/Pq. 

There is also a generalization of the analytic invariant charge Ullfl . Thus, the 
authors of Ref. ^3 proposed the IR enhanced model for the QCD running coupling 



a 



(1) 



*-*-SPPW 



(A2) 



47r 



{z ~ 1) zP 



{a + z - ie)ia + + zP) 



da 



(35) 



where z — j IS? and < p < 1. This formula coincides with Eq. (|ll|l when p = 1. 

Several phenomenological models for the QCD running coupling have been pro- 
posed in Ref. inni It shoul d be rn entioned that the ideas, similar to that of the an- 
alytic perturbation theory^EHI were also used in analysis of the electron-positron 
annihilation into hadronJ^ and investigation of the inclusive r lepton decayl^ 
The thorough study of the power corrections to the strong running coupling was 
performed in Refs. 1721 and E51 There is also a number of methods of the RG im- 
provement of perturbative series for the QCD observables (see, e.g., R,ef. I74|l . 



3. Higher Loop Levels 

3.1. Analytic invariant charge at the higher loop levels 

As it has been shown in Section [21 the one-loop analytic invariant charge possesses 
a number of appealing features. The absence of unphysical singularities and incor- 
poration of the ultraviolet asymptotic freedom with the infrared enhancement in 
a single expression are the most remarkable ones. Apparently, the question of a 
primary interest here is whether these results are affected by the higher loop cor- 
rections or not. So, let us proceed to the study of the analytic running coupling at 
the higher loop levels. 

In accordance with the model proposed (see Subsection 12. l|l . the QCD analytic 
invariant chargeain (g^) at the £-loop level is the solution to the renormalization 
group equation^ 



dln[5ii^(/i^)] 
d In fi^ 



3=0 



J2B,[4'\^.')f\ , 5, = -!^, (36) 



where ai^\fi'^) is the ^-loop perturbative running coupling, a(^^) = Q(/^^)/9o/(47r), 
and /3j are the (3 function expansion coefficients (see [Appendix A| ). By employing 
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the relation the right-hand side of Eq. H36() can be represented as a spectral 
integral 



d In /i^ 



da. 



(37) 



where 



1 



lim VS, 



i=o 



(-CT + ie) 



(38) 



The explicit forms of TZ^^\a) at different loop levels are given in [Appendix C| 
Integrating Eq. H37|l with respect to in finite terms, we obtain 



a 



(?') =«i^„'('Zo) exp 



\l + a/zj a 



(39) 



where z — / and zq = q^/ 

Expression (|39|l contains the explicit dependence on the normalization point q^ 
and on the value of the strong running coupling at this point. It is worth mentioning 
here that this equation possesses the obvious renormalization invariance. In fact, 
having the value of the QCD running coupling at some momenta transferred, one 
might exploit Eq. (|39|) . However, usually it turns out to be more convenient to 
deal with the explicit expression for the strong running coupling independent of 
the normalization point. 

For the case of the model in hand this objective can be achieved by invoking 
the physical requiremenf^ o.al{q^) ai\q^), when ^ oo (see Refs. EHlandlMI 
for details). In fact, this condition has already been employed at the one-loop level, 
see Eqs. (|10ll and Hll() . It is worth emphasizing here that the imposition of such 
requirement does not affect the renormalization invariance of Eq. H39() . In general, 
one is able to derive the integral representation for the QCD analytic invariant 
charge, independent of the normalization point, by invoking the similar requirement 
ain (g^) ain (q^), when ^ oo. For this purpose let us consider the ratio of the 
£-loop analytic running coupling (p??^ normalized at a point q^ to the one-loop AIC 
written in the form H39|l and normalized at the same point q^: 



(^) 



ain^(g2) 



"in (go) 



exp 



\l + a/zfi) a 



(40) 



where Ml^'^\a) = 7^(^)(c^) " 'R-''^\<j). Proceeding to the limit 
at the expression for the ^-loop analytic invariant chargj^^l^ 



oo, one arrives 



An z - 1 



Po zlnz 



cxp 



oc 

jATZ^^^a) ln(l + 



da 



A2- 



(41) 



"This condition has also been used in Ref. 1471 for evaluating the normalization coefficients. 
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Fig. 3. The QCD analytic invariant charge 1411 at different loop levels, z = j . 

It is worthwhile to note that the integral representation of the Kallen-Lehniann 
type also holds for the AIC lUl) 



4^ ^ p^'\a) 



a + z 



da. 



(42) 



where p^^^ (cr) is the f-loop spectral density 
pW(CT) exp 



,0 



dc 

c 



TT 



the one-loop spectral density p^"'^^(cr) has been specified in Eq. H13|l . and 



((t) — TT /A7^W(C) 



c ■ 



(43) 



(44) 



In the exponent of Eq. (j43(l the principle value of the integral is assumed. How- 
ever, in practical use the expression (|41|l turns out to be more convenient. One 
should emphasize here that, similarly to the one- loop case considered above, the 
analytic invariant charge (|41|l has no unphysical singularities, it contains no free pa- 
rameters, and it incorporates the ultraviolet asymptotic freedom with the infrared 
enhancement in an unified manner. 

Figure 13 shows the analytic invariant charge iflTl) ain(g^) — ain(g^) /3o/(47r) at 
the one-, two-, and three-loop levels. It is obvious that AIC possesses the higher 
loop stability. Thus, the curves corresponding to the two- and three- loop levels are 
practically indistinguishable. Proceeding from this one can also draw the conclusion 
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Fig. 4. The /3 function corresponding to the two-loop analytic invariant charge (solid curve). The 
relevant perturbative result is shown as the dot-dashed curve. 

concerning the scheme stabihty of the current approach. In particular, in Fig. |51 
the curve corresponding to the three-loop analytic running coupling a^aniq^) is 
plotted by making use of the coefficient P2 = 2857/2 - 5033nf/18 + 325nf/54 
computed in the MS scheme.I^ The account of the third term on the right-hand 
side of Eq. 10 does not lead to a valuable quantitative variation of its solution in 
comparison with the two-loop approximation. Therefore it is evident that using the 
coefficient (32 computed in another subtraction scheme" does not lead to significant 
variation of the solution to Eq. Q in comparison with the considered case of the 
MS scheme. This statement follows also from the fact that the contribution of every 
subsequent term on the right-hand side of Eq. © is substantially suppressed by 
the contributions of the preceding ones (see also Ref.lMI and [Appendix C|). 



3.2. Properties of the AIC at higher loop levels 

At the one-loop level there is an explicit expression for the analytic invariant 
charge This allows one to perform the investigation of its properties man- 

ifestly. However, at higher loop levels only the integral representation H41|l has 
been obtained for ain('?^) so far. This fact significantly complicates the issue in 
hand, and eventually leads to the necessity of applying the numerical computation. 
Nevertheless, the study of the relevant (3 function enables one to elucidate some 
important questions here, in particular, the asymptotic behavior of the analytic 
invariant charge (see Ref. iMjfor details). 

The /3 functions corresponding to the analytic invariant charge H41|) at the higher 

°At least, in schemes that do not have unnaturally large expansion coefficients (see paper 1761 and 
references therein for the detailed discussion of this issue). 
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Fig. 5. The j3 function corresponding to the three-loop analytic invariant charge (solid curve). 
The relevant perturbative result is shown as the dot-dashed curve. 



loop levels 

d In 



mia)-^-^^^, a(M^) = 5il)(^^) (45) 



are presented in Figures 0] and |S1 Figure 01 shows the two-loop /? function /3an (a) 
together with respective perturbative result (3'^^\a) = —a — Bia^. Figure |31 depicts 
the analogous functions at the three-loop level, namely /3an (a) and /3s (a) = —a — 
Bio? — i?2a^. It is clear from Figures [5101 and [3 that the (3 function corresponding 
to the analytic running coupling H41() coincides with its perturbative analog in the 
region of small values of the invariant charge at any loop level. In other words, 
in the framework of the model under consideration the complete recovering of the 
perturbative limit in the ultraviolet domain takes place. 

Let us turn now to the asymptotics of the /3 function (|45|l at large values of the 
running coupling. FigureElpresents the /3 functions /3in (a), i = 1,2,3 and the one- 
loop perturbative result f3i^\a) — —a. This figure clearly shows the perturbative 
limit at small a, as well as the universal asymptotic behavior of the /3 function 
in hand: pin (a) — —1, when a — s- oo. The latter statement can also be proved 
in an independent way. Indeed, due to the infrared enhancement of the analytic 
running coupling, the value of the right-hand side of Eq. HH6|I (it is nothing but 
the (3 function at the relevant loop level) when — s- corresponds to the limit 
aan(/i'^) — > OO- One can show that irrespective of the loop level 



lim 



{5W(^^)} =1, hm ([5^(^2)1^^') =0, (46) 

I. Jan p'^O {I J Jan 



where j is a natural number (j = 1,2,3, ...). Therefore, we infer that at any loop 
level the /3 function (|45|l corresponding to the analytic invariant charge (|41f) tends 
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Fig. 6. The /3 function corresponding to the analytic invariant charge t41i at different loop levels 
(solid curves). The one-loop perturbative result is shown as the dot-dashed curve. 



to the universal limit 



lim C)(a) = -1. 



(47) 



As it has been mentioned above, such a behavior of the f3 function leads to the IR 
enhancement of the running coupling, namely ~ A^/g^, when — > 0. 

Thus, the /3 function corresponding to the analytic invariant charge H41f) has 
the universal asymptotic behaviors at both small {pin (a) ~ —a, when a —>■ 0) and 
large (/^in (a) — —1, when a oo) values of the running coupling irrespective of 
the loop level. Therefore, the analytic invariant charge itself possesses the universal 
asymptotics both in the ultraviolet (ain(g^) ~ 1/ ln((j'^/A^), when oo) and 

infrared (ainC?^) — A^/g^, when q^ 0) regions at any loop level. In particular, 
this implies that using the AIC (|41|l at the higher loop levels will also result in the 
confining quark-antiquark potential (see Subsection 14. 



3.3. Extension to the timelike region 

In the previous sections the model for the QCD analytic invariant charge has been 
studied in the spacelike (Euclidean) region. However, for congruous description of 
a number of strong interaction processes (for example, inclusive r lepton decay 
or electron-positron annihilation into hadrons) one has to use the continuation of 
the invariant charge to the timelike (Minkowskian) region. In general, one might 
expect such a continuation to affect the properties of the strong running coupling. 
Therefore, it is of a certain importance to extend the analytic invariant charge to 
the timelike domain, and to compare the functions obtained. 

The consistent description of hadron dynamics in spacelike and timelike domains 
remains the goal of many studies for quite a long time. Thus, the first attempts 
to build up the expression for the strong running coupling valid for the timelike 
momenta transferred were made in the early 1980s (see Refs. I77H79|I . Recently, in 
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Fig. 7. The integration contour in Eq. I48i . The physical cut of the invariant charge o(— ^) is 
shown along the positive semiaxis of real f . 

Ref. 1391 the procedure of continuation of the invariant charge from the spacehke to 
the timehke region (and vice versa) was elaborated by making use of the dispersion 
relation for the Adler D function. In the framework of this method the integral 
relationships between the QCD running coupling in Euclidean and Minkowskian 
domains were established. In particular, the continuation of the strong running 
coupling to the timelike region is connected with the "spacelike" invariant charge 
through the relation (see paper ^| and references therein) 

s—ie 

S(.)-^ / ^ = V>0. (48) 

s+ie 

In this equation the integration contour goes from the point s + ie to the point 
s — is and lies in the region of analyticity of the function a{—C). Here and further 
the running coupling in the spacelike region is denoted by a{q^), and in the timelike 
region by a(s). It is worth mentioning also that the relation inverse to (|48|l is 

oc 

J [s + q-'Y 


It is convenient to choose the integration contour in Eq. H48|l in the form pre- 
sented in Figure [7| Namely, the integration path starts from the point s + ie, pro- 
ceeds in a parallel way with the real axis to infinity, then along the circle of infinitely 
large radius it goes counter-clockwise to the point {po — ie), and then it goes in a 
parallel way with the real axis to the point s — ie. Since the QCD running coupling 
possesses the property of asymptotic freedom (namely, a{q^) ~ 1/ ln((7^/A^), when 
q^ —> oo), the integral along the circle of an infinitely large radius gives vanishing 
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Fig. 8. The one-loop analytic running coupling in the spacelike {q^ > 0) and timoliko {s = > 
0) regions, z = q'^ /A'^. 



contribution to the right-hand side of Eq. (|48() . Then, the remaining expression can 



be reduced t 



39. 



ais) = / g{C) 



dC 

C ' 



(50) 



where g{() denotes the spectral density of the KaUcn-Lehmann representation for 
the running couphng a{q^): 



ZTTl e-*0+ I 



(51) 



Thus, the continuation of the analytic invariant charge H41|l to the timelike region 
is given by 



oo 



(52) 



s/A2 



with the spectral density p^^-'(C) denoted in Eqs. H43|l and (|44|l (see Ref. I4(il for 
details). 

The one-loop running coupling H52|l has the following asymptotic in the ultra- 
violet limit: 



N 47r 1 f 1 

"'^^ ^ ~ /3o Inz \ 3 (lnz)2 



O 



1 



1 



(lnz)4' z 



(53) 



where z = s/h?. On the one hand, this running coupling has correct UV behavior 
determined by the asymptotic freedom. On the other hand, the so-called Tr^-terms 
have also appeared in expansion (|53|l . These terms play a key role in the description 
of the strong interaction processes in the timelikc domain (see discussion of this issue 
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in Refs. I38L 1391 1801 and l81|l . It is worthwhile to mention also that, similarly to the 
"spacelike" case, there is the IR enhancement of the one- loop running coupling (|52() : 

^-(^)-ll(l^' ^''^ 

However, here the type of the infrared singularity differs from that of the invariant 
charge Nevertheless, it is this behavior of the running coupling H54(l that 

enables one to handle the integrals over the infrared domain of the form (see, e.g., 
Ref.El 

S 

A{s)^-^Ja,^i.s')ds'. (55) 



It is interesting to note also that within the current approach a simple relation 
holds between the (3 function corresponding to the running coupling in the timelike 
region H52|l and the relevant spectral density: 

d[ail^(s)] 47r 
a Ins po 

Thus, the j3 function in hand proves to be proportional to the spectral den- 
sity p^^\s). Obviously, this result justifies the attempts, originated by SchwingerjSSI 
to find a direct physical interpretation of the (3 function (the so-called Schwinger's 
hypothesis, see also Refs. 15^ and 1511. 

The plots of the functions a^an{q^) and Sin (s) are shown in Figure |H1 In the 
ultraviolet limit these expressions have identical behavior determined by the asymp- 
totic freedom. However, there is asymmetry between them in the intermediate- and 
low-energy regions. The ratio of the one-loop analytic running coupling in the space- 
like region ain ((f') to its continuation to the timelike region ain {—<f') is presented 
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in FigureEl The relative difference between these functions is about several percent 
at the scale of the Z boson mass, and increases when approaching the infrared do- 
main. Apparently, this circumstance must be taken into account when one handles 
the experimental data. In particular, this will play a crucial role in description of the 
inclusive r lepton decay and the e^e~ annihilation into hadrons in the framework 
of the approach developed (see Subsections 14.31 and 14.41 respectively). 

4. Phenomenological Applications 

As has been mentioned above, there are several ways of incorporating the analyt- 
icity requirement into the rcnormalization group method. Eventually, this leads to 
different models for the running coupling in the framework of the analytic approach 
to QCD. Certainly, every such model, as well as any new model for the strong in- 
teraction, must be verified thoroughly. A decisive test of a model's self-consistency 
is its applicability to description of diverse QCD processes. Since we are working 
within a nonpcrturbative approach, the study of intrinsically nonpcrturbative phe- 
nomena is of a primary importance here. Furthermore, the elaborated model for the 
analytic invariant charge contains no free parameters, i.e., similarly to perturbative 
approach, Aq^d remains the basic characterizing parameter of the theory. Hence, 
the congruity of the estimations of this parameter, obtained for different strong 
interaction processes, would imply the quantitative consistency of the developed 
model (see also papers ITTI 1551 and references therein). 

4.1. Static quark antiquark potential 

Let us proceed to construction of the static quark-antiquark potential in the frame- 
work of the one-gluon exchange model. This potential is related to the strong run- 
ning coupling a(g^) by the three-dimensional Fourier transformation 



(see, e.g., reviews 151 II II and references therein for details). In fact, the lowest-lying 
bound states of heavy quark systems may well be described by the potential (|57|l 
even with the perturbative expression for the QCD running coupling However, 
at large distances, which play the crucial role in the hadron spectroscopy, the per- 
turbative approach blows up due to unphysical singularities (such as the Landau 
pole) of the strong running coupling. 

For the construction of the quark-antiquark (QQ) potential we shall use here the 
analytic invariant charge (see also Refs. 1871 1451 1881 andl89|l. After integration 
over the angular variables, Eq. (|57(l acquires the form 



oo 




(57) 







oo 




(58) 
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where p — qro, R = r/ro, and vq is a reference scale of the dimension of length, 
which will be specified below. This equation is rather complicated to be integrated 
explicitly. Nevertheless, one is able to study the asymptotics of the quark-antiquark 
potential H58(l by making use of the following prescription. First of all, let us intro- 
duce a dimensionless variable Q — pR and rewrite Eq. (|58|l as follows: 



16 1 1 

VM = 

^ ^ 3/3o ro R 



R^ 



sin Q dQ 



InE-lnQ 



(59) 



Then, formally expanding the denominator of the integrand, p one can find the 
asymptotic behavior of V{r) at both small and large distances: 



Vir) 
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sin Q 



i?^(lni?)"+iy Q 



(lnQ)"dO 



OO 

mo ^ r ' D 



(60) 



The values of no and nio will be specified below. 

For evaluation of the expansion coefficients in Eq. ()6U|I it is worth considering 
an auxiliary integral of the form 



rfi + -^^ 

Q* sinQdQ = 7^2* ^ ^' 







(61) 



Differentiating this equation n times with respect to variable t we obtain 

f (inQY' sinQdQ = v{n,t), (62) 







where 



r 1 



(63) 



In Eq. (|()T|l the parameter t takes the values < |Re {t +1)1 < 1 (see, e.g., 
R.ef. I90|l . However, in order to determine the expansion coefficients of the second 
sum in Eq. H60|l . one has to go to the point t — —3, which is outside of this range. 
Nevertheless, it can be done by making use of the analytic continuation of the left- 
hand side of Eq. (|61(l . This continuation is unique and it is defined obviously by 
the right-hand side of Eq. 1)61(1 all over the complex i-plane except for the points 



PA similar method has also been used in Ref. |S] 
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t — ~2N, with N being a natural number. Fortunately, we are not dealing with 
these values of the parameter t, and we can put 



sinQ 



(In Q)"dQ = v{n,t) 



, s^2N (iV = 1,2,3,...). 



(64) 



It should be noted here that this analytical continuation plays the role of regular- 
ization of the expansion coefficients in Eq. H6Ufl . It is also convenient to introduce 
the notations 



— v{n, t) 



— w(m, t) 



t=-i 



(65) 



t=-3 



The explicit expressions for these coefficients can be found in [Appendix D| 

Thus, the static quark-antiquark potential H60I) can be represented now in the 
following form: 



^T^s 



R ^ (lni?)"+i 



^ (lni?)™+i 

m— ^ ^ 



. (66) 

7-0 



At small distances the potential H66|) possesses the standard behavior, determined 
by the asymptotic freedom 

Stt 1 1 
3A) n) R Ini?' 



V(r) 



0. 



(67) 



At the same time, it proves to be rising at large distances 



= 7775 



R 



3/3o ro 2 Ini?' 



(68) 



implying the confinement of quarks. It is of a particular interest to mention here 
that a similar rising behavior of the QQ potential has been proposecl^ a long time 
ago proceeding from the phenomenological assumptions. 

As it has been shown in Subsection l3.2l the analytic invariant charge (|41|l has the 
universal asymptotics both in the ultraviolet and infrared regions at any loop level. 
Furthermore, the approach in hand possesses a good scheme stability. Therefore, 
neither higher loop corrections, nor scheme dependence, can affect qualitatively 
the results obtained. The use of different subtraction schemes (or different loop 
approximations) will result in redefinition of the scale parameter, as it also takes 
place within the perturbative approach. 

Equation describes the behavior of the static quark-antiquark poten- 

tial l|58|) at small and large distances. However, its straightforward extrapolation 
to all distances encounters poles of different orders at the point i? = 1, which ap- 
parently is an artifact of the expansion For the practical purposes it would be 
undoubtedly useful to derive an explicit expression for the QQ potential applicable 
for < r < CX3. 
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Fig. 10. The quark— antiquark potential Iti9l in dimcnsionless units at different levels of approx- 
imation: no = mo = (dashed curve), no = mo = 4 (solid curve), no = mo = 5 (a), and 
no = mo = 10 (■)■ The one-loop perturbative result is shown by dotted curve. 



In order to construct an explicit interpolating formula for the quark-antiquark 
potential (|58|l we shall employ here the following method. Let us modify the expan- 
sion 1)66(1 in a "minimal" way, by adding terms which only subtract the singularities 
at the point B. = 1 and do not contribute to the derived asymptotics. This leads to 
the following expression for the static quark-antiquark potential 




R {In RY 



mo 

Reg m=0 
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Reg} 


_(lni?)™+i 





(69) 



where R = r/ro. The functions (lni?)")]Rcg and [i?/(lni?)'"]Rcg are presented 

in an explicit form in [Appendix D| 

The numerical analysis of Eq. I)6H|I revealed that for practical purposes it is 
enough to retain the first five expansion terms {uq = mo = 4) therein. It turns 
out that the potential ((69|l itself and the corresponding estimation of the value of 
parameter Aq^d are not affected by higher-order contributions. In particular, the 
curves ((69|l for no — mo — 5 and for no ~ mo — 10 are practically indistinguishable 
of the curve corresponding to no — mo — 4 over the whole region < r < cx3 (see 
Figure nnj. And the higher-order estimations of the parameter Aqcd vary within 
0.5 % of the value obtained at the no = too — 4. 

Thus, we arriv e a t the following explicit expression for the static quark- 
antiquark potential!^ 
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Fig. 11. Comparison of the quark-antiauark potential V{r) defined by Eq. 1701 (solid curve) 
with the quenched lattice simulation data^''' (•). The values of the parameters are: A = 670 MeV, 
rtf = 0, Vb = — 3.164 GeV. The dashed curve corresponds to the relevant one-loop perturbative 
result. 
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(70) 



where R — r/rQ. In order to reproduce the correct short-distance behavior of the 
QQ potential, the dimensional parameter rg in this equation has to be identified 
with Ams by the relation r^^ = Aexp(7), where 7 ~ 0.57721... denotes the Euler's 
constant (see, e.g., Refs. 1571 1^ and IT^ . It is straightforward to verify that the 
potential (|70|l satisfies also the concavity condition 

^^>0, ^^<0. (71) 

dr dr'^ 

which is a general property of the gauge theories (see Ref. '911 for details). 

Compariso n of the quark-antiquark potential (|70|l with the quenched latt ice 
simulation datcP^ has been performed with the use of the least square method.l^ 
The varied parameters in Eq. H7()|l were and the additive self-energy constant Vq . 
The result of this fit is presented in Figure ^2 This figure shows that in the non- 
perturbative physically-relevant range 0.3 fm < r < 1.2 fm, in which the average 
quark separations ^ (r^) for quarkonia sitsf^SI the QQ potential (|70|l reproduces 
the lattice datJ^ fairly well. At the same time, in the region r < 0.05 fm the de- 
rived potential coincides with the perturbative result (dashed curve in Fig. Illfl . 
The difference between the lattice data and the expression (|70f) in the interme- 
diate range 0.05 fm < r < 0.3 fm can be explained by the presence of additional 
nonperturbative contributions at small distances (see, e.g., Ref. I92|l. 
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Thus, the estimation of parameter Aqcd in the course of this comparison gives 
^(nt=o) _ j-gyg _|_ g-j ivleV (this value corresponds to the one-loop level with rif — 
active quarks) . It is worth noting that the fit with the use of the maximum likelihood 
method also gives a similar value A^'^J"^"^ ~ 650 MeV. Further, in order to collate the 
obtained result with forthcoming estimations, one has to continue it to the region 
of three active quarks. It was performed by employing the matching procedure, and 
givesi the valuc"^ A("f=3) ^ (590 ± 10) MeV. 



4.2. Gluon condensate 

Let us proceed to another nonperturbative aspect of the strong interaction, namely, 
to the gluon condensate. The latter is determined as the matrix element 

K = lim ( vac 



(vac|^:G;;„(a:)G^,(y):|vac), (72) 



the averaging being performed over the "true physical" vacuum state |vac) (see 
Rcf. for a recent review of this issue). In t he f ramework of perturbative ap- 
proach the quantity H72|l identically equals to zero However , lattice simulationJ^ 
and comparison of the QCD sum rules with the experimental data on the low- 
energy hadron dynamic J ^"^ ^ ^ testify to the nonzero value of gluon condensate. 
The quantity (|72|l is the subject of thorough studies for a long time (see, e.g., re- 
views i93l|96j and Refs. 1521 1^1 . In particular, the integral relation between the gluon 
condensate (|72|l and the nonperturbative part of the strong running coupling was 
established in Refs. [TSl 1551 and lUHl 

oo 

J a„p(fc2)fc2dfc2, (73) 



where a^p{k^) — a{k'^) — Q;g(/c^). Here a{k'^) is the QCD invariant charge, and 
ag(fc^) denotes the perturbative running coupling. 

Thus, in the framework of the current approach, the gluon condensate (|73f) at 
the one-loop level takes the form (see Eq. 

oc 



PoTT^ J In z ' 



where z = fc^/A^. However, the integral on the right-hand side of this equation 
diverges. In order to regularize it we employ here the method proposed in Ref. 1961 
Namely, let us introduce the cutoff parameter feg, which has the physical meaning 
of the characteristic scale of the nonperturbative effects. The phenomenological 



"^It is interesting to note here that similar values of the parameter Aqcd have also been obtained 
within different approaches to this issue (see Refs. l5H7t . 
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estimation of this parameter gived^ = (700 ± 100) MeV. Then, assuming the 
integration in Eq. H74|l as the principle value one, we obtain 

where Li(a:) denotes the logarithm-integral^ 

X 

Li(a;)=y^, Rea;>l. (76) 



The solution to Eq. (|75|) with the phenomenological value of the gluon 
condensattPZI K = (0.36 ± 0.02 GeV)" gives A = (631 ± 79) MeV (this estima- 
tion corresponds to the one- loop level with rif = 3 active quarks). This value of the 
parameter Aqcd agrees fairly well with its estimation obtained in Subsection 14. II 

4.3. Inclusive r lepton decay 

In the previous subsections the analytic invariant charge has been applied to de- 
scription of some intrinsically nonperturbative issues of QCD, the congruous esti- 
mation of the parameter Aqcd being obtained. For further verification of the self- 
consistency of the model proposed, it is worth proceeding to the hadron processes 
which are commonly described within the perturbative approach. Investigation of 
such processes at relatively small energies is of a primary interest here. Among 
them, the inclusive r lepton decay is a most sensitive to the infrared behavior of 
the strong running coupling. Thus, we turn now to this hadron process, restricting 
ourselves at this stage to the one- loop level (see Ref. 0Slfor details). 

The experimentally measurable quantity here is the inclusive semileptonic 
branching ratio 

r(r' ^ hadrons" Ur) . 
= _ — T — ■ (77) 

1 (r —fe VeVr) 

One can split this ratio into three parts, namely Rr = Rt.v + Rt,a + Rt,s- The 
terms i?7-,v and Rr.A account the contributions to Eq. H77|l of the decay modes with 
the light quarks only, and they correspond to the vector (V) and axial-vector (A) 
quark currents, respectively. The accuracy of the experimental measurement of 
these terms is several times higher than the accuracy of the strange width ratio i?r,s , 
which accounts the s quark contribution to Eq. 1)77(1 . Thus, let us proceed with the 
nonstrange part of the Rr ratio associated with the vector quark currents (see 
R,efs. .98. .99. .80. and .40. for detailed description of this issue): 

Rr.v = ^ iKdP^EW (1 + <5qcd) . (78) 

In this equation A^c = 3 is the number of quark colors, \V^ ni \ = 0.9734 ± 0.0008 de- 
notes the Cabibbo-Kobayashi-Maskawa matrix element ,1221 = 1.0194±0.0040 
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is the electroweak factorp^^^ and (5qcd is the strong correction. A r ecent exper- 
imental measurement of the ratio ()78(l by ALEPH Collaboration gav^^ -Rt.v = 
1.775 ±0.017. 

In the framework of the approach in hand the one-loop QCD correction is de- 
termined by the integra^^^ 

2 

where (s) is the analytic running coupling in the timelike region l|52f) . Mt- = 
(1776.99l:o:26) MeV deno tes the r lepton mass, and m = (4.0 ± 1.5) MeV stands for 
the light quark mas^^^ at a normalization scale of 2 GeV. It is worth noting here 
that there is no need to involve the contour integration in Eq. (|79|) . since analytic 
running coupling (|52() contains no unphysical singularities in the region s > 0. In 
other words, the integration in Eq. H79|l can be performed in a straightforward way. 
Introducing the notations x — s/M"^, xq ~ Am'^ /M^, and cq = xq{2 — 2xq + Xq), 
one can rewrite Eq. H79|l in a more convenient form 



c _ 1 

Oqcd — 

TT 



aW(M,2) - coSW(4m^)] + ± J (2 - 2x' + x') [x^^ dx, (80) 



where the spectral density p^^^ (g ) is d efined in Eq. 

For the value R^.y given abovJ^^the estimation A("f=2) ^ (55^ ± 69) MeV has 
been obtained for two active quarks. The uncertainty here is due to the errors in 
the values of i?T-,v, iKidl, ■S'ewj itt-, and Mr- However, in order to collate this result 
with the previous estimations of the parameter A opp , one has to continue it to the 
region of three active quarks (the s quark mas J '''' I nis = (117.5 ± 37.5) MeV has 
been employed here). This gives the value A^"''^'^) = (524±68) MeV, which perfectly 
agrees with all the estimations of Aqcd in the framework of the approach in hand. 



4.4. e+e annihilation into hadrons 

It is worthwhile to consider one more strong interaction process, namely, the 
electron-positron annihilation into hadrons. The experimentally measurable quan- 
tity here is the ratio of two cross-sections 

, cr(e+e^ hadrons) 

Unlike the inclusive scmilcptonic branching ratio l|77|l . this equation contains the 
explicit dependence on the center-of-mass energy of the process. In the framework 
of the approach in hand, at the one-loop level the ratio R{s) can be represented in 



32 



the following form (see, e.g., Refs. ll05l andHT |l : 



ni 



R{s) ^N.^Qf 



TT 



(82) 

/=i 

In this equation iVc = 3 is the number of colors, Qj denotes the charge of the /-th 
quark, and Sin (s) is the analytic running coupling in the timelike region H52|l . 
Recent experiment al rn easurement of the R{s) by CLEO Collaboration at 
= 10.52 GeV gavJ^ni Ri^so) = 3.56 ± 0.01 (stat.). With this value the esti- 
mation of the parameter Aqcd for Uf — A active quarks is A*^"'^'*^ = (393 ±65) MeV. 
It is worth noting that this estimation turns out to be rather rough due to a consid- 
erable systematic errors of the experimental data at small energies. Nevertheless, 
the continuation of the obtained result to the three active quarks region gives 
jY(nt=3) — (4QQ ± 76) MeV, that agrees fairly well with all the previous estimations 
of the parameter Aqcd. 

Thus, the applications of the model developed to description of diverse strong in- 
teraction processes ultimately lead to congruous estimation of the parameter Aqcd- 
Namely, at the one-loop level with nf = 3 active quarks A — (557 ± 36) MeV. Ap- 
parently, this testifies that the analytic invariant charge substantially incorporates, 
in a consistent way, both perturbative and intrinsically nonperturbative aspects of 
Quantum Chromodynamics. 



5. Conclusions 

Let us summarize the basic properties of the model for the QCD analytic invariant 
charge, considered in this review, and the main results obtained in this field. 

The proposed way of involving the analyticity requirement into the renormal- 
ization group method results in qualitatively new features of the strong running 
coupling. Namely, analytic invariant charge has no unphysical singularities, and it 
incorporates the ultraviolet asymptotic freedom with infrared enhancement in a 
single expression. Furthermore, additional parameters are not introduced into the 
theory. The consistency of the model developed with the general definition of the 
QCD invariant charge is proved. 

The detailed investigation of the one-loop analytic running coupling and the rel- 
evant f} function is performed. The latter is proved to coincide with its perturbative 
analog at small values of invariant charge. At the same time, this (3 function pos- 
sesses the asymptotic behavior corresponding to the infrared enhancement of the 
strong running coupling. The renormalization invariance of the analytic running 
coupling is shown explicitly. 

The analytic invariant charge is derived and studied at the higher loop levels 
also. It is proved that the higher loop corrections do not affect the AIC qualitatively. 
In particular, the analytic running coupling has the universal behavior both in 
ultraviolet and infrared regions at any loop level, and it possesses good higher loop 
and scheme stability. 
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The analytic invariant charge is extended to the timelike region. It is shown to 
have asymmetrical behavior in the intermediate- and low-energy domains of the 
spacelike and timelike regions. This difference turns out to be rather substantial and 
should be taken into account when one handles the experimental data. The obtained 
result confirms hypothesis due to Schwinger concerning the relation between the 
/? function and the relevant spectral density. 

The developed model is applied to description of diverse strong interaction pro- 
cesses both of perturbative and intrinsically nonperturbative nature. The derived 
static quark-antiquark potential is proved to be confining at large distances. At 
the same time, it has the standard perturbative behavior at small distances. The 
analytic invariant charge reproduces explicitly the conformal inversion symmetry 
related to the size distribution of instantons. The developed model is also applied 
to description of gluon condensate, inclusive r Icpton decay and electron-positron 
annihilation into hadrons. The congruity of the estimated values of the parameter 
Aqcd testifies to the self-consistency of the approach developed. 
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Appendix A. Perturbative QCD Running Coupling 

The QCD invariant charge g^{ix^) is the solution to the renormalization group 
equation 



d In/i^ 



/3(5(m'))- 



(A.1) 



In the framework of perturbative approach, assuming the running coupling being 
sufficiently small, one can approximate the /? function on the right-hand side of this 
equation by the power series 
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(A.2) 
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where for the SU(3) gauge group /?o = 11 — 2nf/3, /3i = 102 — 38nf/3, j32 = 
2857/2 — 5033nf/18 + 32nf/54, and rif denotes the number of active quarks. The co- 
efficients /3o and (3i are scheme independent, whil e the value given for the three-loop 
coefficient /32 is calculated in the MS schemel^ Let us rewrite RG equation HA.1|I 
at the ^-loop level as follows: 



(3, 



Pi 



(A.3) 



where a^in"^) = g^(/x^)/(47r) and a(^2) = a(^2)/3o/(47r). 
At the one-loop level Eq. HA.3|) reads as 



dln[5(^^(^2)] 



d In/i^ 

Integrating this equation in finite terms, one gets 
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In 



(A.4) 



(A.5) 



It is convenient to introduce here the parameter A of the dimension of mass, which 
absorbs all the dependence on the normalization point q^: 



A2 



9o exp 



An 



1 



f^Oai'\ql)\ 



(A.6) 



In this case, the solution to the RG equation (|A.3ll takes the well-known form'' 



/3o In z ' 



A2- 



(A.7) 



Let us note here that this expression for the one-loop perturbative running coupling 
has the unphysical singularity (the so-called Landau pole) at the point q^ = A^. 
At the two-loop level Eq. (jA.3p acquires the form 



d\n[ai'\p^)] 
d In/i^ 



(A.8) 



Integrating this equation likewise the previous case, one arrives at the transcenden- 
tal relation,*' which determines the two-loop perturbative running coupling 



1 



ai^\q^ 



Bi In 



1 



1 



Bia^s'\q^) 



(A.9) 



''It is worth mentioning here that the relation IA.6t ensures the renormaUzation invariance of the 
solution to the RG e quation IA.7I . 

''Certainly, Eq. lA.Si has also solutions different from Eq. IA.9I . Nevertheless, all the ambiguity 
here can be eliminated by redefinition of the parameter A on the right-hand side of Eq. iA.9i . 
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The parameter A in this formula differs from the one-loop one (|A.6 



^(2) = M^exp^ 



An 



/3o aP{^I^) 



Bi In 



1 47r 



Si/3oa(2)(^2) 



(A.IO) 



The explicit solutio n to Eg. (|A.9|1 can be written down in terms of the so-called 
Lambert W functiorP^ (see [Appendix B| | : 



47r 1 



/3o Bi 1 + Wfe{-exp[-(l + Bi-Mnz)]}' 



(A.ll) 



where z = / A? and Bi = (3i/0^. Here the branch index of the Lambert W 
function may take the values k — Q ox k = —1. But only the latter one satisfies the 
asymptotic freedom condition. It is worthwhile to note that two-loop perturbative 
invariant charge ljA.ll(l also has unphysical singularities, namely, the pole at — A? 
and the cut < < A^. The exact solution ljA.ll|l is widely used along with the 
iterative solution to Eq. (|A.9(I 



An 



1 



f3o Inz + Bi ln(l + BfMnz) 
and its approximate form (see, e.g., Ref. ^J): 
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(A.12) 



(A.13) 



Since the account of the higher-order corrections leads to substantial technical 
complications, we present only the explicit expression for the three-loop running 
coupHng here (see, e.g., Ref. llUOfI : 



An 1 



In z 



In z 



Bl 
(lnz)2 



(In(lnz))' 



ln(lnz) + ||-l 



(A.14) 



The calculation of the four- loop (3 function coefficient (3^ in the MS scheme has been 
performed in Ref. ll(J6l For the relevant four- loop perturbative running coupling see, 
e.g., Ref.rnTTl 



Appendix B. The Lambert W Function 

For the representation of the analytic running coupling (|ll|l in the renorminvariant 
form and for the investigation of the corresponding [3 function (I22II it proves to be 
convenient to use the so-called Lambert W function. As long ago as the middle of 
18th century this function is being employed in diverse physical problems .^---- ^ In 
contemporary QCD studies the interest to this function arose just a several years 
ago. Namely, it was revealecP^ that the explicit solution to the perturbative RG 
equation for the invariant charge at the two-loop level can be expressed in terms of 
the Lambert W function (see also [Appendix A| ). 
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Fig. 12. The function xe^ (dot-dashed curve) and two real branches of the Lambert W func- 
tion IB. II : 14^0(3;) (solid curve) and W—i{x) (dashed curve). 



The Lambert W function is defined as a many-valued function Wk{x), which 
satisfies the equation 



Wkix) exp Wfc(x) 



(B.l) 



where k denotes the branch index of this function. Only two real branches of the 
Lambert W function, the principal branch Wo{x) and the branch W-i{x) (see 
Fig. 112(1 . will be used in our consideration. The other branches of this function take 
an imaginary values. One can show that for the branches Wo{x) and W-i{x) the 
following expansions hold: 



Wo{e) 



0, 



W-i{~e) = lne + C'(ln|lne| 
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= -1 



/2^ + C'(e) 



-1 - V2^ + 0(e), 
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(B.2) 
(B.3) 

(B.4) 
(B.5) 



where e = 2.71828 . . . denotes the base of the natural logarithm. Details concerning 
the properties of the Lambert W function can be found in Ref. 11081 

In order to represent the analytic running coupling (|ll|l in a renorminvariant 
form and to derive the corresponding p function ((22(1 one has to solve the equation 
(see Subsection E31 and Refs. l^l^ 

^—^=a, z>0, a>0 (B.6) 
z In z 

for the variable z. Let us multiply this equation through by the factor a^^ Inz 



In z = b, 

z 



z > 0, 



(B.7) 
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where b = —l/a. Obviously, the equation obtained has a trivial solution z — I, 
which does not satisfy the initial Eq. ljB.6|) when a ^ 1, that is a consequence 
of the multiplication of the latter by Inz. Therefore, when solving Eq. (|B.7p for 
the variable z one has to discard its trivial solution z = 1. Next, let us represent 
Eq. HB.7|1 in the form 

^exp(-^=be\ (B.8) 



z \ z 



Taking into account the definition HB.1() . solution to Eq. HB.8|I can be expressed in 
terms of the Lambert W function 

- = Wk{be''). (B.9) 

z ^ ' 

The branch index k of the W function will be specified below. 

In the physically relevant range a > the argument of the Lambert W func- 
tion in Eq. (|B.9|I takes the values — 1/e < b^ < 0. The only two real branches, 
the principle branch Wo{x) and the branch W-i{x) (see Fig. El correspond to 
this interval. Here one has to handle carefully with the interchange between these 
branches at the point x = — 1/e (this corresponds to the value a = 1). Thus, the 
nontrivial solution to Eq. ljB.7|l for the variable z is 



1 ( b-^W^i(be'') , -l<fo<0, 

- = <^ (B.IO) 

^ [b-^Wo{be'') , b<-l 

(another choice of branches will be considered below). Therefore, the solution to 
Eq. ljB.6|l we are interested in can be written in the form 



1 



7V(a)' 

where the function N{a) (see Fig. I13|l is defined by 

r No{a), < a < 1, 
Nia)^{ Nkia)^-aWk 
tiV_i(a), Ka, 



[ 1 ( 




exp 




a V 





(B.ll) 



(B.12) 



As follows from the definition ljB.l|l . the function xe^ and the Lambert W 
function are mutually inverse. Here one has to distinguish precisely the branches 
Wq and W-i (see Fig. I12|) . From this figure it follows directly that for a; > —1 the 
function xe^ "corresponds" to the branch Wq, and for x < —1 the function xe^ 
"corresponds" to the branch W-i. Let us introduce for convenience the notations 
of the inverse functions for these branches 

iWoy\x), x>-i, 

xe''={ (B.13) 
{W-i)-\x), x<^l. 
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0.0 0.5 1.0 1.5 2.0 



Fig. 13. The function N{a) (boldface curve), and the functions No{a) (solid curve) and N-i{a) 
(dashed curve) (see Eq. IB. 121 ). 



Now the function N(a) in Eq. (|B.12|I can be written in a more compact form 



Nia) = 



-aWo[{W-iy\-a-^) 
^aW^i\{Wor'{~a-^) 



< a < 1, 

1 < a. 



(B.14) 



It is worth noting also that another choice of the branches of the function N{a) 



N{a) 



iV_i(a), 0<a<l, 



,iVo(a), 



1 < a 



(B.15) 



leads to the trivial solution of Eq. ljB.7|) (in this case N{a) = 1). Therefore, the 
solution z = 1/N{a) does not satisfy Eq. (jB.6|) when a ^ 1. 

Let us outline briefly the basic properties of the function N(a) introduced in 
Eq. (jB.12|l . In the physical range a > it is a non- negative, monotonously increasing 
function (see Fig. I13() . With the help of the series (|B.2|) - (jB.5|l one can show that 
for the function N{a) the following expansions hold: 



N{a 0+) = exp 
N{1 

N{a ~> cx)) = alna <{ 1 + O 
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(B.18) 
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Appendix C. Spectral Functions 

The spectral function T?.^^^ (ct) at the i-\oop level has the following form (see equa- 
tion ijSHl)): 



7^W(^) = ^ lini y A^J 5W(-a-*e) 



a^^^(-(T + ie) 



(C.l) 



Here ai^'^ {y?) is the £-loop perturbative running coupling and (3j denotes the /3 func- 
tion expansion coefficient (see [Appendix A) ). It turns out to be convenient to rear- 
range equation ljC.l(l : 

n^'\a)^Y.Bjyli{^), (C.2) 

where 



a^^)(-o- + ie) 



(C.3) 



and Bj = l3j/f]l'^^. It is worthwhile to quote also several commonly used relations: 
lim ln(-cr ± ie) = Incr ± ITT, tr > (C.4) 

£—►0-1- 

and 



\n{a±ib) ^ ^\n{a^+b^)±i 



IT fa 

arctan — 

2 \b 



6 > 0. 



(C.5) 



Here it is assumed that — 7r/2 < arctan(a;) < 7r/2. 

The one- loop perturbative running coupling is determined by Eq. ljA.7|) . In this 
case the spectral function (|C.ip takes a simple form 

For the construction of the two-loop spectral function we shall use the expres- 
sion ljA.12|) .'^ which gives 

where 



A(^)(y) 



2 ' 



(C.7) 
(C.8) 
(C.9) 



*It has been shown in Ref. l38l that the spectral functions constructed by making use of Eqs. IA.12i 
and iA.14i practically lead to the same result as the spectral functions corresponding to the exact 
(numerical) solutions to the RG equation lA.Si at the higher-loop levels. 
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with the two-loop functions A*-^^(y) and C^^\y) defined by 



A(2)(y) = l + ^iiln 



C(2)(y) = f + Bi 



y l2 
"f 



1 (y_ 



-\nBi 



arctan 

2 TT V 



Bi 



(C.fO) 
(C.ff) 



Here and further we use the notation y = ln(cr/A^). 

By making use of the perturbative running coupHng (jA.14|l . the foUowing result 
was obtained for the three- loop spectral function l|C.l|l : 



where 



C(3)(y)(3[A(3)(y) 



(y2 +7r2)9 

and the three- loop functions A'''^''(y) and C^'^\y) are 

A(^)(y) - 2/(2/' - 3^2)i?(y) - tt'Ctt' - 3y^)J{y) 

C(3)(y) = 2/(3^' - y')Jiy) - (3y' - ^')i?(y), 

i?(2/) = 2/' - tt' - Bi(2/a - tt'c) + ^'(a^ 

^(y) 



2/ 

2y 
1 



2 2 
TT C 



B, 



Bi_iyc + a) + Bfc{2a-l) 



ln(2/2+7r2), c: 



■ arctan 



(f) 



(C.12) 

(C.13) 
(C.14) 

(C.15) 

(C.16) 
(C.17) 
(C.18) 
(C.19) 

(C.20) 



1 

2 -v. ■ "" " " 2 TT 

The plots of the spectral functions HC.6I) . IIC.71) . and (|C.12ll are presented in 
Figure ITU It is interesting to note that TZ'^^^a) —> TZ^^^a), when ct — > oo, as well 
as when cr ^ 0, irrespective of the loop level. It is worth mentioning also that the 



relation max | fjj (ct) | 3> max | gj^^-^ (cr) | holds at any loop level. In particular, this 
implies that the contribution of the every subsequent term on the right-hand side 
of Eq. (|36() is substantially suppressed by the contributions of the preceding ones. 
In turn, this results in a good higher loop and scheme stability of the approach in 
hand. 



Appendix D. The Expansion Coefficients of the 
Quark Antiquark Potential 

In this appendix the explicit expressions for the expansion coefficients of the static 
quark-antiquark potential generated by the analytic invariant charge H69|l are gath- 
ered in a concise form. By making use of Eqs. H65|l and (|63|l we find 
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-15 -10 



10 15 



Fig. 14. The spectral functions H^'^^ct) defined by Eq. icTTl . The one-, two-, and three-loop 
levels are shown by the dot-dashed, solid, and dashed curves, respectively. 
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where 7 0.57721... denotes the Euler's constant, and ((x) is the Riemann's Zeta 
functiorP^ (C(3) — 1.20206...). Applying the "regularization" procedure specified 
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in Subsection 14. II one obtains (see Eqs. H69() and (|70|l 'l 
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